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ABSTRACT 


An experimental investigation of the internal ος 
of solid fuel ramjets was conducted in order to determine 
the regression rate of the fuel as a function of chamber 
ши-ввиге, inlet air temperature, and air flux rate, and to 
model the flow in solid fuel ramjets which use sudden expan- 
sion flameholders at the inlet. In addition, flame stabili- 
zation nios wepeerdnvestisated. A computer solution for 
the non-reacting flow field gave results in good agreement 
with experiments. 

solid fuel pramjets have an average regression rate of 
the fuel that closely follows the theoretical expression 
derived for kinetically controlled hybrid rocket combustion. 
The inlet step-height-to-motor-diameter ratio is the domi- 


nant parameter in determining flame stability limits. 


However, inlet velocity is also an important parameter. 
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IS S ENUPRODUCTION 


The solid fuel ramjet can be classified as a special 
case of the hybrid rocket. That is, the fuel is solid and 
the oxidizer is gaseous (air). The solid fuel ramjet pro- 
vides many of the advantages of solid fuel rockets without 
the weight penalty of carrying its own oxidizer. Ease of 
handling, storability, weight savings, and possibly greatly 
increased tactical ranges are among its advantages. 

A research project is currently being conducted at the 
Naval Postgraduate School to model the internal ballistics 
of solid fuel ramjets. Such a model would provide the 
design engineer with a working tool that would permit 
Mmiesormance parameters and operating characteristics to be 
Predicted. The overall goal of the research project is to 
velop a model that can predict the regression rate of the 
fuel as a function of operating and configuration variables. 
E сте experimental data currently exist that can be 
used to determine the adequacy of a model, thus the model 
must be developed in shone verifying each step with experi- 
mental data.  Non-burning tests to determine flow patterns, 
Don burning tests with wall blowing to determine the effect 
of wall mass Een the flow characteristics, and hot 
meine tests with various inlet geometries are required. 
This thesis initiates the overall research project 


objectives. 





Information is provided from several areas of study 
that are pertinent to the internal ballistics of solid fuel 
Bets. Specific areas are flow and heat transfer in. ducts 
uua backward facing step, flame stabilization in high 
speed flows, hybrid rocket combustion, and thermal decompo- 
sition of polymers. 

Abbott and Kline [Ref. 1] investigated the turbulent 
flow over backward facing steps using water flow through a 
e i mensional channel. It was found that the recircula- 
pusnczone consisted of a complex pattern involving three 
distinet regions. Immediately downstream of the step face 
there is a three-dimensional zone "characterized by one or 
Mire vortexes rotating about an axis normal to the vertical 
wall." This zone is followed by an almost two-dimensional 
zone of recirculating flow. The area at the flow reattach- 
ment point is characterized by a fluctuating three- 
dimensional flow. It was also found that the location of 
the reattachment point remained fixed as Reynolds number 
and/or inlet een intensity varied, and that increas- 
ing the step height moved the reattachment point further 
downstream from the step face. | 

Krall and Sparrow [Ref. 2] using water flow through an 
orifice in an electrically heated eircular tube found that 
local heat transfer coefficients in the separated, reattached, 


and redevelopment regions are several times larger than those 


for fully developed turbulent pipe flow. They found that 





SE maximum heat transfer coefficient occurred in the 
mily of the point of flow reattachment. 

At low Reynolds numbers based on downstream pipe. 
diameter (10,000 to 50,000) the heat transfer near the 
reattachment point varied rapidly with Reynolds number. 
However, at higher Reynolds numbers (70,000 to 130,000) 
ENS weak dependence was noted. They also found that for 
weaker . (small step heights, h) the reattachment 
point was well defined but for strong separations (large 
step heights) the reattachment point was spread out. This 
transition from a single well-defined reattachment point to 
the "spread-out" configuration occurred between h/D=0.16 
and h/D=0.25, where D is the downstream pipe diameter. The 
reattachment point was located between 1.25 and 2.5 pipe 
diameters downstream of the initial separation. Increased 
TES height moved the reattachment point slightly downstream. 
Reuewds number and Prandtl number had no significant effect 
on the reattachment point location. 

The results of Krall and Sparrow are also in general 
agreement with the results of Abbott and Kline. For small 
step heights (h/D=.17) the two-dimensional and axisymmetric 
reattachment points are very nearly equal (X,/D from 1.070 
1.25). However, for larger step heights, the reattachment 
zone spreads out considerably more for the two-dimensional 
case than for the axisymmetric case. 

The ideal hybrid rocket is one in which a gaseous 


oxidizer flows over a solid fuel. In a hybrid rocket, gas 
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phase combustion in the form of a diffusion flame takes 
place in the turbulent boundary layer formed by the oxidizer 
flow over the solid surface of the fuel. Heat is trans- 
ented by convection and radiation to the solid surface 
which causes decomposition of the fuel. The convective heat 
Io cusport is complicated by the presence of high rates of 
Wall mass addition. Using a heat transfer limited model, 
regression rates of the fuel are generally proportional to 
E. Ere 5-15 спе total gas Flux and n isma constant 

Goi cally 0.8). 

Experimental data obtained at low chamber pressure 
and/or high oxidizer mass flux indicate that the regression 
IN hybrid grains becomes less dependent on G, but is 
directly related to the partial pressure of the oxidizer in 
the gas phase [Refs. 3 and 4]. Reduced reaction rates at 
lower pressures broaden the flame zone. Schlieren studies 
have also shown that the flame zone broadens at low chamber 
pressures and may broaden enough to undergo transition to 
esturbulent КОСО С: flame [Ref. 6]. When the oxidizer is 
diluted, the diffusion flame moves further away from the 
іссе, increases іп thickness, and the flame temperature 
drops [Refs. 3 and 5]. In addition, surface kinetics may 
be rate controlling under these operating conditions 
[Ref. 7]. When these conditions exist, the heat transfer 
limited model does not adequately describe the combustion 


process. 


Tk 





Hybrid rockets operating at high chamber pressures are 
adequately described by a heat transfer limited model. 
Muzzy [Ref. !] has reviewed the modifications that have been 
mS lo the heat transfer limited model to account for gas 
phase and surface kinetic effects. The modified models 
[Refs. 6 and 7] predict a dependence of the regression rate 
on pressure as well as mass flow rate. For a system where 
radiative heat transfer can be neglected, Muzzy has shown 

\ 


The regression rate equation in the kinetically domi- 


L region of operation can be expressed as: 


το 


where C is a constant, P is pressure, G is total mass flux, 
eS The longitudinal coordinate, a is a constant that 
antes from 0.3 to 0.4, and b is a constant that varies from 
-0.1 to -0.2. Another study by Kustov and Rybanin [Ref. 18] 
νι. τ: treated the pressure sensitive regime of operation. 
In a solid fuel ramjet the reacting flow has character- 
istics that are similar to those discussed above. However, 
two significant differences exist. Solid fuel ramjets would 
normally operate at lower chamber pressures than hybrid 
rockets, resulting in significantly greater port velocities. 
In addition, the oxidizer in this ease is air which would 
correspond to a severely diluted oxidizer in hybrid combus- 
tion. These differences require some means of flame stabi- 
lization. A series of tests performed at the United 


Technology Center [Ref. 8] have shown that one means of 
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sustaining combustion in a solid fuel ramjet is to use a 
sudden expansion inlet as is used in dump combustors. 
Similar results were reported by McCarthy [Ref. 9]. 

NWucchychoveloetty gas stream, where the velocities of 
the flow are at least an order of magnitude greater than the 
turbulent flame speeds of the combustion mixtures, flame 
stabilization is required. Bluff body flame stabilization 
k cen used successfully in turbojet and liquid fuel 

' 
ramjet combustor design in the past. Typical flame stabi- 
lizers have been V-gutters and rods, or the dome of a can- 
type combustor [Refs. 10 and 11]. In addition, dump or 
sudden expansion combustors have been studied for integral 
rocket/liquid ramjet designs [Ref. 19]. 

The recirculation zone behind the bluff body or step 
ο σος intense mixing of fresh unburned fuel and air with 
hu combustion products. The hot combustion products ignite 
Nel alr mixture. This stabilized flame can then 
τε throughout the combustor. Due to the intense 
PreoulLent mixing Та ο ο ο ο πα оп Zone.  COMDUSTION may 
in some сасес approach that of a wetrt-stirred reactor. 

A well-stirred reactor model considers the composition 
ana all of the thermodynamic properties to be uniform 
throughout the volume. An idealized reactor with instanta- 
neous mixing would have homogeneous fuel, air, and product 
concentrations, and constant pressure and temperature 


Mets. 12 and 13]. 
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Solid fuel ramjet combustors have, depending upon the 
combustor length and inlet configuration, two separate and 
distinct combustion zones, each of which appears to be 
controlled by a different mechanism [Ref. 9] (See Figure 1). 
In the recirculation zone products and reactants are rapidly 
med, and in the limiting case a well-stirred reactor may 
be approached. Downstream of the flow reattachment point 
the boundary layer develops and the combustion may be similar 
CoO that | a hybrid rocket in which kinetic effects are 
dominant (low pressure, high mass flux, severe oxidizer 
ше топ). 

Thus, one possible model for solid fuel ramjet combus- 
tion would be to treat the recirculation region as a well- 
stirred reactor and the downstream region as a kinetically 
dominated hybrid. The larger the grain length-to-diameter 
ratio, the more dominant the downstream region becomes in 
determining the average regression rate, and the more one 
expect the ramjet to behave as a kinetically controlled 
hybrid ШО оздор. However, for small L/D ratios and/or 
large Beeireulation zones, the combustion in the recircula- 
tion region may become important in determination of the 
КЕ герге regression rate. In addition, the combustion char- 
wei Stics in the recirculation region determine the 
flammability limits for the entire combustor. Thus, the 
model used in this region requires considerable attention. 
fer very strong recirculation Zones, a well-stirred reactor 


model may be adequate. However, for weaker recirculation 
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Memes, Which are desirable from a pressure loss standpoint, 
kinetic rates may be of the same order or greater than the 
reactant mixing rates. In this ease, a well-stirred reactor 
model cannot be expected to be applicable and a more detailed 
Esuncs4deration of the recirculating flow is necessary. Thus, 
the work of Gosman, et al. [Ref. 14] and Spalding, et al. 
{Ref. 15] becomes increasingly important and may provide 
the basis for modeling the recirculation region in solid 
fuel ramjet combustors. 

It is the purpose of this study to determine the 
dependence of the regression rate of a solid fuel ramjet 
we hamber pressure, inlet air temperature, and mass flux 
EE. and to model the flow in the recirculation region 


caused by the sudden expansion inlet. 
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II. METHOD OF INVESTIGATION 


Non=burning flow visualization tests were conducted 
mE a Plexiglass model with three inlet step configura- 
mmn. flow reattachment points of the incoming stream 
were measured for various step Reynolds number flows. 

A Fortran language computer program, Pistep II, devel- 
oped by Spalding, et al. [Refs. 14 and 15] provided a finite 
difference solution of the two dimensional Navier-Stokes 
equations governing heat and mass transfer in recirculating 
flows. It was used to calculate the flow reattachment point 
and the flow field for non-burning flows. The results were 
compared with the flow visualization data and the work of 
previous investigators. The boundary conditions on the 
differential equations were altered to simulate blowing at 
the wall in order to study the effect of mass injection 
Шеш спе wall оп (пе flow field. 

Experimental firings were conducted using polymethyl- 
methacrylate grains Dee with non-vitiated air to provide 
regression rate data as a function of inlet air temperature, 
chamber pressure, and oxidizer flux rate. These tests were 
made to provide initial data for characterizing ramjet 
combustion and flame blow-off limits, and to determine if 
the reacting flow behavior was similar to that obtained 


analytically and experimentally for non-burning flow. 
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ПО  DESCRIFTIOT OF APPARATUS 


A. COLD FLOW VISUALIZATION MODEL 

The flow visualization model (Figure 2) consisted of 
three sections. The first section was a 2% x 2% x 5-inch 
plexiglass block with a three-quarter inch diameter hole 
bored through the center. It received air from high pressure 
tanks and supported a perforated disk flow straightener. 
The Step section was made of stainless steel and was ten 
inches in length to further dampen turbulence and provide 
lE developed velocity profiles. Three steps were used 
with inside diameters of 0.246, 0.519, and 0.714 inches 
respectfully. The third section was a 15%-inch plexiglass 
tube with an inside diameter of 1.469 inches. 

mer the lower flow rates a sonic choke with a diameter 
of 0.081 inches was used to measure flow rates. For higher 
flow rates, an orifice flow meter with inputs to a Visicorder 
was used to measure flow rates. The flow reattachment 


points were measured visually using a scale. 


B. HOT FIRING APPARATUS | 
A schematic of the solid fuel ramjet apparatus showing 
the overall system is presented in Figure 3. 
1. Ramjet Motor 
The solid fuel ramjet motor consisted of four main 
sections: the head end assembly, the step insert section, 


the grain, and the aft closure section (Figure 4). A11 of 
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the metal parts of the motor except for the bolts were made 
ENEStalnless steel. 

The head end assembly was 6% inches long and contained 
τ ἵπ1ε6ετ5 Tor the air, the oxygen-methane ignition mixture, 
E mitrogen purge, and the cooling air. Two Champion Z7 
EE plugs were used for ignition. 

The step insert section was 2-3/4 inches long and 
held the step in place. Three steps were used during this 
meso ol cation with inside diameters of 0.25, 0.50, and 0.75 
inches. This section also accomodated the forward end of 
the PMM grain. 

τος Ansthacnylate was chosen as the fuel in 
this study because it has been widely used for basic research 
in hybrid rocket combustion and in studies of polymer 
er eradation in the past. It is composed of long chains of 


Sealing monomer: 


H E 

| | 
P unm 

| | 

H COOCH 


Idealized chemical reaction with complete combustion 


would involve a stoichiometric reaction with air as follows: 


C HO, + 60, + 22.56N, > HHO - 2005 T 22. 50N, 


> 


The average molecular weight of the products inthis case 
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is 31.3 lbm/lbmole. The gas constant is R/M = 49.4 
ft-lbf/lbm-?R, and the stoichiometric air fuel ratio is 

8.27 lbmair/lbmfuel. The equilibrium ratio of specific heats 
at a temperature of 3500°R is y = 1.24. An adiabatic flame 
memeecracvure calculation for the stoichiometric reaction using 
a value of Ho ӨТ 92.8 keal/mole for the heat of formation of 
PMM at 25 °C [Ref. 5] gives a flame temperature (with no 
dissociation) of 4265 °R. 

The polymethylmethacrylate grains were cut from 
a slab into 3-3/4 x 3-3/4 x 12 inch blocks. Тһе ends were 
rounded to an outside diameter of 3-1/2 inches to fit into 
Step Section and the aft closure section, and were 
sealed with "O" ring seals. The initial inside port diameter 
Р ЕПС ргаіпс наз 1.5 inches. Eight grains were made that 
Here six inches in length, but the majority of the tests 
were performed using grains that were 12 inches in length. 

The aft closure section contained a straight 
section with a pressure tap and a converging-diverging 
nozzle. Nine oc les with various throat areas were used 
to control the chamber pressure. A Wiancko 0-200 psig pres- 
sure transducer was mounted on the forward part of the aft 
closure section to measure the chamber pressure. 

The grains were held in place by four 3/8 inch 
diameter steel rods between the step section and the aft 
closure section. Mounting flanges on the head end assembly 
and the aft closure section were bolted to a large test 


stand. 
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2. Ignition System 


The ipnition system metered oxygen and methane at 
800 psig from two high pressure tanks through needle valves 
into the head end assembly section. Inlet ports on the top 
ші πε bottom injected methane, and two inlet ports 
injected oxygen from either side. Two inches downstream 
from the oxygen-methane ports were two Champion Z7 spark 
plugs located 180 degrees apart. A transformer supplied 
power to a Model T spark coil to energize the spark plugs. 
С ivrogen Purge and Cooling Air System 
mme mitrocen purge and cooling air system was 
Em elted upstream of the head assembly section. Two high 
pressure tanks fed nitrogen into the motor at the end of a 
B ring for approximately one second to extinguish the 
Eembusulon flame. Cooling air from a low pressure compressor 
was then blown through the motor. 
ПЕ Air Flow Control 
A standard ASME orifice flow meter was used to 
measure the flow C of the air into the motor. A manually 
operated ο ο ας ρου σε ος Orifice and the motor was 
used to provide the desired flow rate through the motor. 
Two pneumatically operated Jamesbury ball valves operating 
together either vented the air to the atmosphere, or allowed 
me alr to pass through the motor. 
5. Data Acquisition System 
A Colvin 0-35 psi differential pressure transducer 


was used to measure the pressure drop across the orifice 
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flow meter. А Wiancko 0-200 psig pressure transducer was 
uo measure the pressure upstream of the orifice, and 

an iron-constantan thermocouple was used to measure the 
temperature at the orifice. A second iron-constantan thermo- 
couple was used to measure the temperature of the air enter- 
Ing the motor. A Wiancko O-200 psig pressure transducer 
measured chamber pressure at the aft closure section of the 
MOCOT. 

The orifice thermocouple output was recorded on a 
0-600°F strip chart recorder. All of the other transducer 
outputs were recorded on a Honeywell Model 2106 Visicorder. 
A five cycle per second timing signal and an ignition pulse 
were also recorded on the Visicorder. 

6. Air Feed System 

A Pennsylvania air compressor provided air at pres- 
ЕШШЕс up to 150 psia at a rate of up to 700 scfm. This air 
was fed into an air reservoir and a Polytherm air heater. 
The Polytherm air heater was capable of yielding 1.75 pounds 
per second airflow at 150 psia and 1000%F. The heater 
burned natural gas, me Che alr passed through a heat 
exchanger, thus providing non-vitiated hot air. 

A temperature controller Кз са two mixing valves 
that produced the required air temperature at 150 psia. A 
Шива! gate valve controlled the air flow to the test cell. 
A valve in the test cell controlled the flow of air to the 
pipe leading to the orifice flow meter. 

Figures 5, 6, 7, and 8 show the test cell and the 


Ramjet apparatus. 
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' ERIT MENTAL PROCEDURE 


BEE COLD FLOW TESTS 

The cold flow tests were performed at the rocket lab 
[15165 οἱ the Naval Postgraduate School. Visualization 
was obtained by injecting a mixture of water, alcohol, and 
iss coloring at the base of the step inlet through a small 
tube with a hypodermic needle. The mixture was broken up 
Into E... small droplets by the turbulence in the recircu- 
lation zone. In the recirculation zone the droplets on the 
wall moved towards the step. At the flow reattachment point 
che droplets on the wall were stationary. Downstream of the 
reattachment poínt the droplets on the wall moved away from 
ши-щетсер. Using this technique, the flow reattachement 
point locations for various step Reynolds number flows using 
three step sizes were measured visually using a scale. 

Initial visualization tests were attempted using a 
tufted model but they proved to be of limited value due to 


Hn rohly turbulent nature of the flow. 


БЕЙНЕСІ FIRING TESTS 

All hot firing tests were performed in the jet engine 
test cell at the Naval Posteraduate School. Hot firing 
data were obtained from a series of firings at three alr 
flow rates, three chamber pressures, and three inlet temper- 
atures, using three step heights and two grain lengths. The 


majority of the hot firings were performed with a 12 inch 
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grain length and an inlet diameter of 0.5 inches. Table I 


summarizes the nominal test conditions. 


TABLE I 


NOMINAL TEST CONDITIONS 


Chamber Pressure (psia) 50 15 100 
Air Flow Rate (lbm/sec) πη 5 .40 
Inlet Temperature (°F) 70 200 350 
Inlet Diameter (in) 25 „50 ЛЕ 
Step h/D „417 Ecc DE 
Grain Length (in) 12 6 


For the tests in which the Polytherm air heater was used, 
the temperature controller was set and the temperature would 
stabilize at the ramjet inlet in one to two hours. This 
delay was due to the fact that the flow rates to the motor 
were quite low, at most 0.5 lbm/sec. Once the temperature 
Peel zed au tne desired temperature at the inlet to the 
ramjet, the manual gate valve downstream of the orifice was 
set to obtain the desired air flow rate. 

Chamber pressures were preselected and controlled by 
Varying the nozzle throat area for a given air flow rate. 
Chamber pressures could normally be attained within ten psi 


of the desired pressure. Calculated combustion temperatures 
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varied significantly with air flow rate and step size. Also, 
reicher chamber pressures the preset flow rate would 
decrease upon ignition because the higher back pressure would 
umnchoke the air flow control valve. 

TeEnilutlon normally lasted for four seconds. After three 
Ends of ignition, the air was directed through the ramjet 
meer. Ignition was continued into the first second of the 
run to insure that combustion would be sustained. A series 
of thirteen tests at high flow rates in which the ignition 
flame was blown out were used to determine the rate of 
eonsumption of the PMM grains during the oxygen-methane 
ignition. These data were used to correct the initial 
αἱ the fuel used in the regression rate calculation. 

Combustion lasted for thirty seconds. At the end of 
each run the motor was extinguished by simultaneously venting 
the air to the atmosphere with the Jamesbury ball valves and 
actuating the nitrogen purge system. Low pressure air was 
mS blown through the motor for cooling. The only section 
of the motor that heated up significantly was the aft 


closure section and the nozzle. 
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E ES EE τ COMPUTER PROGRAM 


Pistep II is a Fortran computer program developed by 
Gosman, Pun, Runchal, Spalding, and Wolfshtein [Ref. 14]. 
It provides a finite difference solution of a steady 
incompressible turbulent flow through a cylindrical pipe 
with a sudden step enlargement. Equations [Ref. 15] 
Ex crpning the distribution of stream function, y, vorticity, 
ш/т, ME imensional temperature, T, turbulent kinetic energy, 


LEE uU turbulent energy dissipation rate, е, are cast in 


ENNIO form in cylindrical polar coordinates as follows: 


CE a) 
9 9 9 7 9 $ 
E A за J 
: 3(C, 9 
- —[b,r ' cec “O 


where q represents any one of the variables v, w/r, k, €, 


C апаи οτε presented if) 


EI The functions a b $? ф 


2 vo 
Table II. 
A two equation turbulence model is employed in the 
ри обра. The two variables are turbulent kinetic energy, k, 
and the turbulent energy dissipation rate, e. The effective 
wescosity, Hoff? is related to k and e by 


E 2 
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A "dissipation length scale", 2, can be used in an 


alternative formulation: 


ша 


Mere = SK PA 


Башшар the two equations for the effective viscosity gives: 
L = x? je 


Constants that appear in the two turbulence equations 


are C4, C5, C, , сү, and σε" Values used in the computer 


u 
program together with their Fortran symbols are given in 


Table III. 
TABLE TI 
VARIABLE FUNCTIONS (Table I of Ref. 15) 
a b C S 
$ $ $ $ $ 
QD ` po a Hoff 0 
| ° 
Y 0 1/pr a -w/r 
k l p 1 -(G-pe) 
2 
€ 1 μο ρρ/σς 1 -(c]Ge/k - сәре /К) 
T 1 μρρρ/ση 1 0 





G is a generation term: 
9 V д ð V 
= Zy2 42 
G = Mantel + 0 за 1! 


V V ον 
ri2 r\ı2 7, 
kuwa <= <= с ее Е 
T ) + {ο DNE bo 
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TABLE TIT 


TURBULENCE CONSTANTS (Table II of Ref. 15) 


Constant C1 Co “| σι. σ. 
Fortran Cl C2 CMU PR(NK) | PR(NE) 
Symbol 

DU 1.45 2.0 0.09 1.0 ἵν; 


Boundary conditions are specified at every point on the 
boundary. Wall functions that are derived from Couette flow 
solutions are used in evaluating boundary values for w/r, k, 
ουσ Other boundary conditions are given in Figure 9. 

The computer program calculates heat and momentum 
transfer in steady two-dimensional axisymmetric recirculating 
flows. Five Simultaneous coupled non-linear partial differen- 
tial Bons are solved. The finite difference equations 
р сіс еапасіолз that are coupled and non-linear. The 
Gauss-Seidel point iterative method is used to solve the 
finite difference equations. Provisions are made for both 
БЕ рге |ахаС1оп апа underrelaxation of the variables. 

References 1! and 15 contain a more complete discussion 
of the mathematical formulation of the problem, development 
of the finite difference equations, and solution algorithms. 
An excellent discussion of several turbulence models can be 


mound in Reference 16. 


e 





The boundary conditions on the pipe wall can be changed 
so that the effect of mass addition through the wall for 
non-reacting flows can be investigated analytically. A 
constant radial velocity at the wall, which yields a 
linearly varying stream function along the wall, is used to 
simulate wall mass injection (See Figure 10). In this 
initial study, the main flow and the wall mass addition 
were identical in chemical composition and temperature. 

in Pistep II, wall functions based on Couette flow are 
used to specify the boundary conditions on the pipe wall 
for all of the variables except stream function. For very 
low blowing rates, which are characteristic of solid fuel 
ramjets, Couette flow solutions may yield realistic results. 
However, for higher mass injection rates, Couette flow 
solutions may not adequately represent the boundary 
conditions. 

The centimeter-gram-second system of units is used in 


ENMCNMPQIstep II computer program. 
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I πεις AND DISCUSSION 


A. COLD FLOW EXPERIMENTS 

A series of cold flow visualization tests were conducted 
to determine the reattachment point location for flow over a 
rearward facing step at the inlet of a cylindrical tube. 


Defining a step Reynolds number as: 


= p Vh where 


Re step Т 


p fluid density 
V... average velocity in the port of the tube 
Е Step height 


ee lua viscosity 


it was found that the reattachment point location remained 
constant as the step Reynolds number was changed for a given 
SCEp height. The reattachment point moved further downstream 
with increasing step height. The tube inside diameter was 
used to non-dimensSionalize the reattachment point location 
from the face of the step. A linear least-square curve fit 
of 47 data points was used to obtain the following equation 


(See Figure 11). 


x 
T = 7.15 (h/D) + 0.35 


Figure 12 compares the data obtained in this investigation 


ГЕР the data of Abbott and Kline [Ref. 1] and Krall and 


3 






Sparrow [Ref. 2]. The unsymmetric reattachment points 
obtained by Abbott and Kline are shown. The flow reattached 
at different locations on either side of the two-dimensional 
channel. The location of the maximum heat transfer coeffi- 
pent found by Krall and Sparrow for water flow through ап 
orifice in a cylindrical tube was used to locate the 
reattachment point data shown in Figure le. 

The data from this study fell between that of Krall and 
sparrow and the downstream reattachment point found by 
Abbott and Kline. Air was used in this study whereas water 
was used in the other studies; however, each study found 
that reattachment point location was not affected by changing 
Reynolds number for a constant step height. Krall and 
Sparrow also found that changing Prandtl number of the flow 
KOL change the reattachment point location. 

The spread in the cold flow data in Fig. 11 is probably 
due to larger relative errors in measurement closer to the 
step face and not due to Reynolds number effects. 

At higher step Reynolds numbers small vortices were 
observed at the step face similar to the observations of 
Abbott and Kline for two-dimensional channel flow. The 
reattachment point was axially symmetric for most Reynolds 
number flows. However, at high Reynolds numbers the reattach- 
ment point varied around the circumference of the tube. Four 
distinct lobes could be seen when this occurred and may be 


associated with the small vortices at the step face. It 
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was also observed that the rotation and mixing of the flow 
EN Le reorrculation zone beeame more intense as the height 


of the step increased. 


ERA? ESTEP II RESULTS 

ihe Pistep II computer program was run for four step 
Configurations and five inlet velocities. Table IV summarizes 
the reattachment point results. A linear least-square curve 
fit is shown in Figure 11 and Figure 12 for comparison with 
L cola flow data and the data from other investigations. 
From these figures it is seen that Pistep II gives a depen- 
dence of ΧΡ on h/D which is in general agreement with 
experimental results. With the current turbulence model and 
wall functions, the results agree closely with those of 
' Krall and Sparrow and fall somewhat below the data obtained 
in this investigation. 

ШЕРІ x 21 grid was used for the h/D of .25 and a grid 
of 20 x 20 was used for the other values of h/D. The 
computing time required for each run was on the order of 20 
emo minutes. Underrelaxation of the variables was required 
to ensure numerical stability. Even with underrelaxation, 
the case with h/D of 0.421 was numerically unstable for 
most Of the inlet velocities tried. However, this case 
Converged the fastest when numerically stable. 

At each of the grid points, Pistep II calculated the 
aAmstribution of the following parameters: vorticity, stream 


function, turbulent kinetic energy, turbulence dissipation 
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* 


xx 


h/D 


.250 
.250 
250 
.250 
. 289 
‚289 
‚289 
„342 
Eu 
„342 
„342 
„421 
421 


Based upon a fixed laminar viscosity of 1.78 x 10 


TABLE IV 


RESULTS OF PISTEP II CALCULATIONS 


X, /D 


55 
1.696 
10001 
ο 
ΠΠ 
2.008 
2.004 
СІБЕ 
2.199 
2161 
2.215 
2.705 


| 2.708 


gm/em.sec 


Re step 


2446 
12232 
24465 
36698 
2009 
20089 
31033 
1335 
6677 
πο ο 
20032 
4109 
6163 


* Niter** V4 , Cem/sec) 
900 1524 
700 7620 
900 15240 
700 22860 
800 1524 
800 15210 
800 22860 
700 1524 
800 1920 
700 15240 
800 22860 
ӘС 15240 
Sal 22860 


H 


Niter is the number of iterations 


за 





meee. mon-dimensional temperature, axial velocity, radial 
r ocity, effective viscosity, density, and turbulence 
length scale. Variations of shear stress and heat transfer 
coefficient along the walls were also calculated. The 
location of the flow reattachment point was determined by 
solving for the point at which the wall shear stress was 


Zero. 


ο PISTEP II WITH WALL MASS ADDITION 

The boundary conditions were altered along the "north" 
wall as shown in Figure 10. A linearly varying stream 
function (i.e., a constant radial velocity at the wall) 
was used to simulate injection of mass uniformly through the 
La τον this initial study, the chemical composition of 
the mass injected through the wall was chosen to be identical 
to that of the main stream coming in the inlet. FAR is a 
term that is analogous to a fuel air ratio, and represents 
the fraction of the total incoming inlet stream that is 
injected through the wall. For example, FAR = 0 gives no 
wall mass prd dou FAR = 0.5 gives an amount of wall mass 
injection equal to half of the inlet stream so that the amount 
leaving at the exit is 1.5 times the nom entering at 
mae inlet. 

FAR was varied in increments of 0.1 for the case in 
which inlet velocity was 500 ft/sec (VIN = 15240 cm/sec) 
foreman П/Г СІ 0,3042, Stream function output was punched 


onto cards which were used as input for a separate program 
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eener rated Contours of streamlines using the Calcomp 
plotter (see Appendix A). Figures 13 through 21 show the 
Streamline patterns obtained for various conditions. It 
should be noted that the streamline patterns are shown only 
for tche upstream half of the ramjet and that the radial 
coordinate has been expanded more than the axial coordinate 
Eu (пе figures аге not to scale. This permits a better 
visual presentation of the streamlines in the recirculation 
zone. Note that with wall mass addition (Figures 14 through 
18), a streamline that originates on the north wall in the 
melemmeulatvion region escapes from the recirculation zone and 
enters the main stream of flow. Figures 20, 13, and 19 show 
the effect of increasing step size (h/D) while maintaining 
inlet velocity constant (mass flow rate decreases). Compari- 
son of Figures 20 and 21 shows the effect of decreasing 
inlet velocity (or inlet mass flow) with a fixed h/D and 
comparison of Figures 13 and 21 shows the effect of decreasing 
step size with a fixed inlet mass flow rate. 

Table V summarizes the results for wall mass addition. 
The reattachment pom location, based on the point where 
wall shear stress is zero, moves upstream slightly as wall 
blowing is increased. The "wall shear stress" is calculated 
from the stream function difference between the wall node 
moore che first node out from the wall. In the tube between 
these two nodes the axial velocity is at least an order of 
magnitude greater than the radial velocity. This permitted 


the calculation of the "reattachment point" in the same 
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manner as for no blowing. Also, the magnitude of the radial 
velocity at the wall is quite small compared to the magnitude 
of the inlet velocity. It should be emphasized that these 
pesults were obtained using the same boundary conditions on 
w/r, T, e, and k as for no blowing. Further studies are 
poured to determine whether or not this assumption is 


realistic. 


TABLE V 


PISTEP II WITH WALL MASS ADDITION (h/D = 0.342) 


Vin (cm/sec) FAR x, /D Vua] (m/sec) 
15240 0.0 πη 0.0 
15240 pon ES nn 
15240 | Do 2.168 9.5 
15240 0.3 2.145 14.2 
15240 | 0.H 155 19.0 
"15240 0.5 201 a 





Note: 900 iterations for each of the above runs 


D. PMM/AIR RAMJET 
Seventy-four hot firing tests were conducted. A Fortran 
language computer program was written to reduce the data. 


A computer listing together with the output is included in 


>> 





ο  ποιχ DB. The first three pages of the output tabulate 
the measured data for each test in which combustion was 
sustained. This is followed by the calculated results. 
This program was modified to reduce the data for the tests 
Inewhieh combustion was not sustained. A computer listing 
ου Ори for the non-sustaining tests is included in 
appendix C. 

The expected uncertainties in the experimental results 
was calculated using the error analysis method of Kline and 
McClintock [Ref. 17]. Using the notation defined in the 
pet data reduction computer program, the uncertainties 
are given in Table VI. 

The step inlet caused non-uniform regression of the fuel 
Pramas. There was less regression of the fuel at the head 
Eun the recirculation region, and nearly uniform regres- 
EIE MEdownstream of the reattachment point, where regression 
Mamas observed to decrease slightly with increasing 
Eu nce from the head end. The inlet effect was more 
pronounced for larger values of h/D. For the majority of 
the tests there was a significant amount of carbon deposited 
ОП The face of the fuel at the head end. The reattachment 
point locations were measured to the edge of the carbon 
deposits or to the location of maximum regression. These 
two locations were nearly coincident in most cases. The 
regression patterns for three tests are shown in Figure 22. 

Regression rates were calculated using two methods. 


Based upon the exit diamter, average regression rate can 
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TABLE VI 


ERROR ANALYSIS 


Symbol % Error Variable 

WDOT 6.3 air weight flow rate 

το 13.83 calculated chamber 
temperature 

RDE 1.4 r based on exit diameter 

RWT 5 r based on weight loss 

WF .6 fuel weight flow rate 

FG 8.6* gross thrust 

FISP 8.6* fuel specific impulse 

RESTEP 5.0 step Reynolds number 

GAIR D air flux rate 

UPORT 15.0* average port velocity 

SFC | 8.6* specific fuel consumption 

EQN 0 empirical r equation 


The uncertainty in the ratio of specific heats, y, and 
the gas constant, R, were not included in this analysis. 
Complete expansion through an ideal nozzle was assumed 
in calculating FG, FISP, and SFC. 


БД 





be defined by: 


ЕЕ ο ο 


where а, апа de are the initial and final exit diameters 
before and after firing, and At is the burn time. An 
average regression rate based on weight loss of the fuel 


can also be calculated: 
\ 


lj (W, = W.) 
ο ж s "M nuseo 
2 At ue 








where W, "Кате CNC INCI and stanallweights of the 


f 
Eu o is the density of the fuel, and L is the grain 
EU This expression for the regression rate gives а 
beer average value when there is non-uniform regression 
along the length of the grain, and was used throughout 

chis investigation. 

Ine dependencêée of the regression rate of the fuel on 
chamber pressure, P (psia), inlet air temperature, T (°R), 
and average mass flux of air, G (1bm/in*sec), was determined 
by holding two of these parameters at fixed values and 
1ης the third TR tests. It was found that 
си ee pression rate varied as: 


0.51 т0.34 0.41 


p = СР G (in/sec) 
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where C is a proportionality constant. For a step h/D of 
ша В, СЕ 2.3 х 107" and for a step h/D of 0.416, 
с = x ο. A plot of the regression rate based on 


weight loss, Y versus this empirical regression rate 


wt? 
equation is shown in Figure 23. On the sixth page of the 
ramjet data reduction computer program output in Appendix B, 


the numerical results are tabulated. The expected uncer- 


tainty in the empirical equation was calculated to be 11%. 


Muzzy [Ref. 4] has shown that the regression rate 
equation for a hybrid rocket operating in the kinetically 


dominated region can be expressed as: 


where a varies from 0.3 to 0.4 and b varies from -0.1 to 


-0.2, and X is the longitudinal coordinate. 

Thus, the average regression rate of the fuel in a solid 
fuel ramjet is very similar to that for pressure sensitive 
hybrid rocket combustion. This is not surprising since the 


reattachment point normally occurred at less than one-fourth 
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oí the grain length. Aft of the reattachment point combus- 
Should approach that of a kinetically controlled hybrid 
combuústor. 

The maximum regression rate of the fuel occurred near the 
Rosstron of the flow reattachment point. This is due to 
το Гассогв. More oxidizer reaches the wall from the 
impingement of the main oxidizer stream at the reattachment 
Em Ти addition, the maximum Nusselt number occurs near 
the reattachment point location. 

Data from 36 of the tests in which combustion was ΠΟΙ 
sustained are tabulated with the computer program in Appen- 
dix C. For the smallest step h/D used (0.25), combustion 
Was Sustained only while the ignition system was operating. 
As soon as ignition was terminated, the flame blew off. 

For the middle step (h/D = 0.333) the flame blew off as 
Soon as the air flow was activated for many of the tests 
at higher flow rates. Some of the tests in which combus- 
las not sustained may have been due to insufficient 
ion prior to activating the air flow. 

Mest s performed at the United Technology Center [Ref. 8] 
.. imaicated that inlet configuration is more important 
in determining when combustion will be sustained. A critical 
step h/D ratio was found to be the best criterion. For 
values of h/D less than the critical value, combustion could 
not be sustained. It was reported that for one particular 
all-hydrocarbon fuel the critical h/D was 0.14. For 


polymethylmethacrylate they found that the critical h/D 
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DN This criterion did not vary with step Reynolds 
number. 

The data from this investigation indicated that for 
high inlet velocities with h/D greater than the critical 
πο determined by UTC, combustion may not be sustained. 
It appears that not only is h/D important, but inlet 
Wemociuy, ipnition time, and inlet temperature may also 
be Ec canc. At identical air mass flow rates, some of 
the tests at higher inlet temperatures failed to sustain 
combustion. This is probably a velocity effect since at 
constant mass flow rates, increasing temperature increases 
the velocity. 

Ehe Pistep II results, although for non-burning cases, 
should indicate trends to be expected in the burning tests. 
The data indicate that the magnitude of the velocity at a 
peace in the recirculation zone is directly proportional 
to the inlet VOU CUL. Also the velocity in the recircu- 
lation zone V, was found to be inversely proportional to 
h/D. Thus V « VIN and V « 1/(h/D). It was also found 
that E ος Ë but did not vary with inlet velocity. 

A simplified model Sr We wi Lamce stabi Pizavion region 
Bone in which the recirculation zone is considered as an 
ignition zone whose length, L, corresponds to the reattach- 
mn point location. At the flame blowoff condition, the 
Bases traveling through the recirculation zone have an 
average contact time equivalent to the ignition time of 


the gases: ο Ξ ZE e a characteristic velocity 
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in the recirculation zone, then t, « L/V. Since 
© 75 

and 
L « X, /D « h/D 


Шеп one possible expression for contact time is 


om 


h 
D 





The step h/D has a more dominant effect in determining 
flammability limits, however inlet velocity is also impor- 
tant. This treatment is oversimplified and was discussed 
ONE To indicate that the inlet velocity may be of signifi- 
mme ih determining blowoff limits. In actuality, ignition 
vk а function of the fuel-air ratio in the recirculation 
ο πα this is in turn a function of h/D. Thus, the func- 
tional form becomes more complex as one considers the 

es of the recirculation region. 

It should also be noted that the simplified model 
described above neglects the effects of turbulence intensity 
lame stabilization. Examination of Pistep II output 
indicates that for a constant inlet mass flow rate, increasing 


step height increases both the turbulent kinetic energy and 
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| E Turbulence dissipation rate. Also for a fixed inlet 
Ev bv. increasing step height decreases the recirculation 
puesveloecity, the turbulent kinetic energy and the turbulence 
dissipation rate. 

ere аге three blow-otf situations to be considered. 
The first one is a steady-state blowoff limit that occurs 
иеп the mass flow rate in the recirculation Zone is greater 
than the mass consumption rate by chemical reaction. The 
mass flow rate through the zone will depend upon the size 
of the zone and the velocity through the zone. This may be 
mportaat in throttling applications during normal combus- 
EE Ihe second Situation is one in which there are igni- 
tion transient blowoff limits during starting. Low chamber 
E ure during ignition gives higher inlet velocities which 
may cause the flame to blowoff even though combustion would 
Sustain if ignition occurred at the expected operating 
pressure. ΤΟ it is possible that the flow at the 
step inlet an Degome schokedmat highimass flow rates when 
the chamber pressure is low during ignition. Subsequent 
КО Sion could modify the size and flow characteristics 
of the recirculation zone and cause the flame to blowoff. 
This indicates that if a ramjet operating with high flow 
RECS could be started at a higher chamber pressure, 
combustion may be sustained. 

ΠΕ TISTE II computer program should prove to be 
advantageous in studying blowoff limits. When chemical 


reaction and wall mass injection are included in the 
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n ans blowofi velocities in the recirculation zone could 
be realistically studied. The effect of h/D, inlet velocity 
(or mass flow rate), and chemical reaction rate could be 
бш егеда. Blowofif limits have been calculated in this 
manner for premixed gases in a dump combustor [Ref. 15]. 

If the results of Abbott and Kline [Ref. 1] are similar 
to what may be expected in an axisymmetric Situation with 
Combustion, inlet turbulence would not be expected to affect 
the recirculation попе o емет πες игй епее шау 
aifect the regression rate of the fuel downstream of the 
reattachment point. Further tests are needed to clarify 
Eu ес of inlet turbulence and distortion on solid 


s ranjet combustion. 
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VII. CONCLUSIONS 


The cold flow visualization data were in general agree- 
ment with previous studies. Reynolds number had no measur- 
able effect on the size of the recirculation zone and the 
reattachment point location. The reattachment point moved 
further downstream from the face of the step as step height 
bor ased. 

The Pistep II computer program appears to predict with 
acceptable accuracy the reattachment point locations for 
Nol reacting flows with no wall mass addition.  Pistep II, 
modified for wall mass addition, indicated that the size of 
LDEEeereclrculation Zone and the reattachment point location 
aremot greatly affected by wall mass addition. The location 
of the reattachment point moved slightly closer to the step 
face as wall blowing was increased. Experiments are required 
to verify this result, and the assumptions concerned with 
the boundary conditions used in the computer program require 
further ЕА. The Pistep II computer program has 
proven to be accurate enough to warrant further development 
Әр а ballistics program for the recirculation region. 

Solid fuel ramjets which use PMM as the fuel have an 
average regression rate that follows the expression: 


MEC Up 99551 4 0.34 5 0.1 xn addition, the maximum 
regression occurred near the reattachment point location. 


Downstream of the reattachment point the regression rate 
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decreased slightly with increasing distance from the head 
end. Thus, solid fuel ramjets have an average regression 
rate that closely follows the theoretical expression derived 
lxx gneutlcally eontrolled hybrid rocket combustion. 

The step h/D is the dominant parameter in determining 
flame stability limits; however, inlet velocity is also an 
mm@pertant parameter. 

The results of this investigation indicate that a 
reasonable internal ballisties model for solid fuel ramjets 
may be obtained by using the Pistep II computer program 
(modified for chemical reactions and wall mass addition) 
w eher with the regression rate equation for kinetically 
''ττο1]εα hybrid rocket combustion. The modified Pistep 
II program should be used to analyze the recirculation 
region (blowoff limits, etc.) and to provide the input 
data (chemical composition, etc.) to the regression rate 
equation aoe kinetically controlled hybrid rocket combus- 
tion. the latter should be employed downstream of the 


reattachment POING: 
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OTT CES TTONS EOR FUTURE WORK 


This work was an initial study of the internal ballis- 
tics of solid fuel ramjets. Additional experimental and 
analytical studies are required in order to a) obtain 
increased combustion efficiency, b) refine the internal 
ballisties model and с) study the effect of other inlet ge- 
ometries'and flow conditions on the regression rate of the 
fuel. Sie experimental studies needed include deter- 
mination of a) the effect of wall mass addition on the 
recirculation zone for a step inlet using inert flows (to 
compare with the theoretical predictions using Pistep II), 
b) the effect of different inlet configurations such as 
swirl. vanes and plug inlets on fuel regression rate and 
flame stability limits, and c) the effect of different 
mixing techniques on combustion efficiency. Analytical 
studies required are a) modifications to the Pistep II 
computer program to include chemical reaction with heat and 
mass transfer at the wall and b) combination of the recir- 
culation zone model (as input) with the kinetically controlled 
hybrid combustion model to effect a usable internal 


ballistics model for solid fuel ramjets. 
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FIGURE 10. BOUNDARY CONDITIONS WITH 
WALL MASS INJECTION 
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CURE 11. REATTACHMENT POINT VS. STEP HEIGHT 
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Krall and Sparrow - based on location of maximum 
heat transfer coefficient (axisymmetric flow) 


Abbott and Kline - note different reattachment poine 
locations for each side of two-dimensional channel 


COLD FLOW RESULTS 


er mer [T RESULTS 


FIGURE 12.  REATTACHMENT POINT DATA COMPARISON 
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